Explicit exact expression for the Thomas precession 



Sergei A. Klioner 

Lohrmann Observatory, 
Dresden Technical University, 
Mommsenstr. 13, 01062 Dresden, Germany 



GAIA-CA-TN-LO-SK-004-2 
issue 2, 14 March 2008 

This work gives an explicit exact expression for the Thomas precession arising 
in the framework of Special Theory of Relativity as the spatial rotation resulting 
from two subsequence Lorentz boosts. The final result for the orthogonal matrix 
of Thomas precession is given by Eqs. (|2ip - (|25p . A trivial calculation leads to the 
compact formula (|26p for the angle of rotation due to Thomas precession. 

In the framework of Gaia the special-relativistic Thomas precession is an important 
step in the derivation of an aberrational formula with the Mansouri-Sexl parameters. 
The latter formula will be used to test the Local Lorentz Invariance with Gaia data 
as will be explained elsewhere. 
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INTRODUCTION 



The Thomas precession naturally arises in Special Theory of Relativity as the additional 
rotation to be added to a Lorentz boost to represent the result of two subsequence Lorentz 
boosts. Although the derivation of the Tho mas prec e ssion fro m the Lorentz transformations 
can be found in many textbooks (see, e.g., Jackson! (|l975[) orlMollerl ( 19721 )). it is normally 
done in the form of expansion in powers on 1/c. ISalingarod (119861 ) gives the exact expres- 
sion for Thomas p recess ion, but in a form that is not readily useful for further calculations. 
Sexl &: Urbantkd (120011 ) have also given the exact expression, but have not simplified it al- 
gebraically, leaving the reader with a rather lengthy calculations. The purpose of this short 
note is to derive the exact and fully simplified expression for the Thomas precession directly 
using two subsequence Lorentz transformations and representing them as a Lorentz trans- 
formation plus a spatial rotation. All calculations have been performed explicitly and in 
normal vector notations. The resulting formula for Thomas precession (Eqs. (I2~T1) — (I25p be- 
low) is rather compact and is valid exactly. Using this expression for the Thomas precession 
a trivial calculation leads to Eq. (1261) for the angle of rotation due to Thomas precession. 

The exact formula for the Thomas precession is interesting by itself, but can also be 
considered as a st ep in the discussion of th e Thomas precession in the framework of Mansouri- 
Sexl test theory (IMansouri fc Sexllll977l ). That latter discussion is important to interpret 
the results of various modern project performing high-accuracy directional measurements 
(e.g., Gravity Probe B, Gaia or SIM) in terms of the Local Lorentz Invariance (see, e.g., 
Klionerl 120071 . for the case of Gaia and SIM). 
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The notations of this paper are usual: c is the velocity of light in vacuum, lowercase 
Latin indices take values 1, 2, and 3 and refer to spatial components of corresponding 
quantities, index is used for time components, Greek indices take values 0, 1, 2 and 3 
and refer to all space-time components of corresponding quantities, repeated indices (both 
Latin and Greek ones) imply Einstein summation rule irrespective of their positions (e.g., 
di hi = ai&i + 0262 + 03^3), the spatial components of a quantity considered as a 3- vector are 
set in boldface (a = a 1 ), the absolute value (Euclidean norm) of a 3- vector a is denoted as 
a or \a\ and is defined by a = \a\ = (a 1 a 1 + a 2 a 2 + a 3 a 3 ) 1 ^ 2 , the scalar product of any two 
3-vectors a and b with respect to the Euclidean metric 8^ is denoted as a ■ b and defined 
as a ■ b = Sijd 1 V = a 1 b\ the Kronecker symbol (unit matrix) is denoted as <5 lJ , parentheses 
surrounding a group of indices denote symmetrization (e.g., = | (A u + A 7 *)), brackets 
surrounding two indices denote antisymmetrization (e.g., = | (A l i — A^ 1 )). 



II. TWO SUBSEQUENT LORENTZ TRANSFORMATIONS AND THE 

THOMAS PRECESSION 

Let us consider three inertial reference systems: (X° = cT,X l ), (x° = ct,x l ) and (x° = 
ci,x l ). The velocity of x a with respect to X a is V\ The coordinates X a and x a are related 
by a Lorentz transformation of the form 

x a = A a pX> 3 (1) 

where 



A° = r, (2) 

A° a = -TK a , (3) 

Al = -TK\ (4) 

r = (i -K-Ky 1/2 , (6) 

K = -V. (7) 

c 



The inverse transformation reads 



X a = A^ xP, (8) 



where A^j is equal to A^ with —K substituted for K (V remains the same after this substi- 
tution). The velocity of reference system x a with respect to x a is v l , and one has 



x 



a = \%x^ (9) 



where has the same form as A^j with 



7 =(l-fc-fcr 1/2 , (10) 

k = -v (11) 

c 
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substituted for T and K, respectively. Now, the velocity of x a relative to X a is V. Using 
standard considerations one gets the relation between the three velocities: 



K 



1 



1 +p 



r k 



i 



r 



i + r 



P K 



where 



C 

p = k ■ K. 

Combining (CQ) and ([9]) one has the relation between x a and X a : 
Now, let us define matrix kg with the same structure as A? but with 



I K K 



-1/2 



T° - A 

nab \ b 

where P ab is the orthogonal matrix describing the Thomas precession 

pa 



A- 

C- 
V- 



6 ab + AK a K b + B k a K b + C K a k b + V k a k b 

(i-7)r 2 
(i + r)(i + f)' 

7 r (, , n T-7T 



l + r 

7 r 



1 + 2 



1 + 7) (i + r) 



i + r' 
7 2 (i-r) 

i + 7)(i + f) 



(12) 

(13) 
(14) 



(15) 
(16) 



(17) 



substituted for T and K for K. According to ffTTI) . ffTU]) . (jH]) and f[T2"|) one gets 

f = 7 r(i+ P ). (is) 

Straightforward calculations show that 



(19) 
(20) 



(21) 
(22) 

(23) 
(24) 
(25) 



Matrix P ab is orthogonal and satisfies the relation P ac P bc = 5 ab . If k is parallel to K (that 
is, for k = a K with an y a), it is easy to check from (12"T|) - fl25j) that the Thomas precession 
vanishes and P ab = 5 ab . ISexl fc Urbantkd (120011 ) have derived this result, but have not given 
it in explicit and fully simplified form. 
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The angle of rotation a due to Thomas precession can be directly computed from the 
trace of matrix P ab using the standard formula 1 + 2 cos a = P aa . Trivial calculation leads 
immediately to 

(i + 7 + r + f) 

l + cosa = — ^ (26) 

(i + 7 )(i + r)(i + r) 

The equivalent r esults hav e been derived after lengthy calculations by iMacfarlane (Il96l 
Eq. (124)) and Urbantke ( 199Cll . the last equation of the paper) and discussed also by 



Sexl fc Urbantkd fl200ll . Eg. (2.10.7)). 



III. IMPORTANT LIMITS FOR THE THOMAS PRECESSION 



From this matrix one can easily restore all standard results concerning the Thomas pre- 
cession. Expanding P ab in terms of k = \k\ one gets 



P 



ab 



5 ab + 



2T 

i + r 



fc[« K b] + 0(k 2 



(27) 



where A^B^ = | (A l B^ — A^B 1 ) is the antisymmetric part of A 1 B^ for any two vectors A 
and B. Defining 5K = K — K and using (fl2l . Eq. fT27|) can be re- written as 



yab 



5 ab + 



2T 2 

i + r 



SK [a K b] + 0{\5K\ 2 ). 



(28) 



This latter form can be found, e.g., in I Jackson! (119751 ) and iMollerl (119721 ). Finally, expanding 
(I2~l~]) - (l25l) in powers of 1/c one gets 

pab = S ab + k [a R b] + k [a R b] ^2 + R 2 _ ^ 



-- (k 2 K a K b + K 2 k a k b - 2 p k {a K b) ) 
+0(c- 6 ). 



(29) 



This expansion can be conveniently used for modelling of high-accuracy directional data. 
Let us note that the symmetric part of the terms of order (9(c~ 4 ) immediately follows from 
the antisymmetric terms of order 0(c~ 2 ). Indeed, considering a general representation P ab = 
fiab _|_ Y2^ =1 e k f£ b with any formal parameter, the condition of orthogonality P ac P bc = 5 ab 

allows one to determine the symmetric term at any order of e. In particular one has 

f (ab) = Q) f (ab) = _i f M f M f <pb) = _1 ^[oc] f M + f M f [a c] y ^ g f [ab] = fc[a Rb] ^ 

immediately restores the symmetric part of the terms of order (9(c~ 4 ) given in the second 
line of ([29]): e 2 f!f b) = -§ (k 2 K a K b + K 2 k a k b -2pk { - a K^) . 
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